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Abstract 

We show that extended self-similarity, a scaling phenomenon firstly observed 
in classical turbulent flows, holds for a two-dimensional metal- insulator tran- 
sition that belongs to the universality class of random Dirac fermions. Devi- 
ations from multifractality, which in turbulence are due to the dominance of 
diffusive processes at small scales, appear in the condensed matter context as 
a large scale, finite size effect related to the imposition of an infra-red cutoff 
in the field theory formulation. We propose a phenomenological interpreta- 
tion of extended self-similarity in the metal-insulator transition within the 
framework of the random /3-model description of multifractal sets. As a nat- 
ural step our discussion is bridged to the analysis of strange attractors, where 
crossovers between multifractal and non-multifractal regimes are found and 
extended self-similarity turns to be verified as well. 
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I. INTRODUCTION 



The surprising observation of extended self-similarity (ESS) in fully developed turbulence 
[1], a hidden scaling behavior which takes place in the crossover between inertial and dissi- 
pative ranges, has provided an empirical - and very efficient - way of determining accurate 
scaling exponents of velocity structure functions. The essential assumption of the method 
is that there are no anomalous corrections to Kolmogorov's 4/5 law, as required from the 
existence of a constant energy flux through the inertial range [2]. To state in an explicit way 
what is meant by ESS, let the average 

S q (r) = (\v(x 1 )-v(x 2 )\ 9 )= lim ^ C dt\v{x u t) - v(x 2 , t)\" (1.1) 

T^oo 1 Jo 

denote the q-th order structure function, where r = \xi — x 2 \ and v(x,t) is the velocity 
field. In the inertial range, the scaling result S q (r) ~ holds. At smaller length scales, 
where dissipation effects come into play, S q (r) loses its power law form. However, a plot of 
G q (r) = \n[S q (r)\ versus G p (r) = \n[S p (r)}, shows in a clear way that the linear relation 

G p (r) = C ^G q (r) + c(p,q) (1.2) 

is reasonably verified even in the crossover towards the dissipative scale. Since Kolmogorov's 
4/5 law predicts £(3) = 1, improved evaluations from experimental data can be obtained for 
C(p), using q = 3 above. Later on, it was realized that a refinement of ESS, referred to as 
"Generalized Extended Self-Similarity" would lead to better linear fits [3,4]. Actually, one 
finds that if G n (r) is taken as fixed "reference structure function" , for some arbitrary n, then 
G q / n {r-) = G q {r) — (q/n)G n (r) satisfies to a higher degree of precision a relation similar to 
(1.2), with slope given now by [n((p) — p((n)]/[n((q) — q((n)]. From now on we will specify, 
only when necessary, the former definition of ESS, by "ESS-I" while its generalized version 
will be named "ESS-II". 

It is clear ESS-I leads to ESS-II, but the converse may be not true. The definition of 
ESS-II is indeed necessary to cope with the fact that at the smallest scales one expects to 
have S q (r) ~ r q . It is not difficult to prove that ESS-II implies structure functions may be 
written as 

S q (r) = c,I/i(r)]«[/ a (r)]««) , (1.3) 

where /i(r) — > 1, f 2 (r) ~ r in the inertial range and fi(r) ~ r, f 2 (r) — > 1 when r — > 0. 
At present, there is no theoretical understanding of the simple crossover behavior (1.3) for 
the structure functions in a turbulent fluid. We note, however, that in the slightly different 
context of the turbulent transport of a passive scalar, an analytical proof of ESS-I was 
carried out in ref. [5]. 

The idea that ESS could hold for a larger class of models is supported not only by 
the problem of the transport of a passive scalar, but also by numerical investigations of 
kinetic roughening [6] and magnetohydrodynamic turbulence [7]. It is tempting to look 
for ESS whenever multifractal scaling laws take place, adventuring beyond the turbulence 
arena where it has been observed. This is precisely the task we pursue in this work for a 
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few additional models. Our motivation is twofold. We stress that the phenomenological 
description of multifractality in turbulence, based on the random /5-model approach to the 
Richardson cascade [8], suggests that ESS could be analogously found in quantum problems 
like condensed-matter localization (see sec. II). Furthermore, we have in mind the conjecture 
that several multifractal systems could be mapped into specific patterns, likely through 
a field theory point of view [9]. Here, inspiration comes from the sucessful applications 
of field theoretical methods in critical phenomena, where universality classes and critical 
properties of statistical systems at thermodynamic equilibrium proved suitable for systematic 
study. In analogy with the renormalization group strategy, firmly based on the experiment, 
we are now faced, thus, with the primary problem of recognizing, among the plurality 
of multifractal models, common features - like ESS - that could be derived from some 
general theory. Of course, an important role in this scenario is to be played by a number of 
well-established multifractal field theory models, as the ones for random spin systems [10], 
Anderson localization [11], the two-dimensional metal-insulator transition [12-14], hadron 
jets in QCD [15], etc. 

While studying some completely diverse models, we have structured the paper to be as 
self-contained as possible. In the next section we discuss a two-dimensional metal-insulator 
transition, where we respectively verify ESS-I and ESS-II for strong and weak disordered 
regimes of the cutoff field theory. We study the transition via its modeling in terms of 
two-dimensional Dirac fermions [12-14]. A "cascade description" of the metal-insulator 
transition at the critical point is addressed, which takes us naturally to the analysis of 
strange attractors, the subject of sec. III. We study then the two perhaps most celebrated 
strange attractors, Henon [16] and Lorenz [17], collecting clear evidence for the existence 
of ESS-II and ESS-I, respectively, in these geometrical sets. In sec. IV, we put forward an 
analysis of ESS in terms of the random /3-model description of multifractality, which, albeit 
essentially phenomenological, throws some light on the previous findings. In sec. V, we 
comment on our results and point directions of further research. 

II. THE TWO-DIMENSIONAL METAL-INSULATOR TRANSITION 

The multifractal nature of delocalized states at the metal-insulator transition in con- 
densed matter systems has been attracting a large deal of interest in recent years [18]. We 
propose here a qualitative cascade picture to approach the physics of delocalized multifractal 
wavef unctions. A more quantitative, although phenomenological, description of the "cascade 
picture" of the metal-insulator transition will be postponed to sec. IV. 

Suppose that at initial time a smooth and localized wavepacket is defined, centered 
around an arbitrary point P. As the system evolves, quantum propagation in the disordered 
background will make the wavepacket to spread and fragment, so that at very large times 
and very far from P, a cascade for the probability density develops, yielding a multifrac- 
tal measure. Of course, multifractality is lost at scales of the order of the system's size, 
where the smoothness of the wavefunction is recovered due to reflections at the boundaries. 
Such a crossover between two scaling regimes is expected to be found in any experiment or 
realistic model of electron localization. The above state of affairs is analogous to the one en- 
countered in turbulence, if the roles of infrared and ultraviolet regions are exchanged. Eddy 
fragmentation in turbulence gives rise to an energy cascade, and dissipation at smaller scales, 
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where viscous effects are more important than convection, render the velocity configurations 
smooth. 

Our prototype for the study of the multifractal-nonmultifractal crossover of probability 
density functions will be an exactly solvable model which describes a two-dimensional metal- 
insulator transition, closely related to the one observed in the integer quantum Hall effect 
[12-14]. The model consists of Dirac fermions moving in a box of size L x L under the 
presence of a transverse random magnetic field. Deviations from multifractality follow from 
the imposition of periodic boundary conditions in the x and y directions, introducing an 
infra-red cutoff in the corresponding field theory model. The random Dirac hamiltonian is 
given by 



M=l,2 



(2.1) 



where is a random gauge field and the a^s are Pauli matrices. 

The exact multifractal spectrum computed by Castillo at al. [14] is obtained from the 
zero modes of (2.1). The normalized solution of Htp = is just 



ipo(x) = M 1 exp[-0(f)o- 3 - ix(x)] 



(2.2) 



where the gauge field 1 has been written as A^(x) = J2^ e ^ud u (l>(x) + <9 M x(^) an d 



J d 2 xexp(—2(f)(x)) 



(2.3) 



The random magnetic field B(x) = — V 2 0(a 7 ) is logarithmically correlated, that is, its fluc- 
tuations are described by the probability density functional 



P[0] oc exp[-i- / rf 2 x(V0) 2 ] . 



(2.4) 



We are interested to study amplitude fluctuations of ipo(x), which do not depend on the 
phase x(x)- Let Q be a square box of size r x r. We define now structure functions of order 
q as 



S q (r) = ([J n \H 2 d 2 x 



(2.5) 



The above average is computed in a quenched disorder scheme from the probability density 
functional (2.4). According to the exact results [14] we have S q (r) ~ r( r( ^ +2 ) for r <C L, 
where r(q) is the scaling exponent given in table 1. Strong and weak regimes are defined, 
respectively, by q c < 1 and q c > 1, where q c = J2n/g. 



Note that the total magnetic flux on the two-dimensional space is zero. Exact results were 
obtained so far only in the trivial topological sector. 
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Structure Function Order 


r(q) (weak disorder) 


r{q) (strong disorder) 


\q\ < 1c 


2(q-l)(l-q/qi) 


-2(1 - q/qcf 


( 


1 > 1c 


2q(l -sgn(g)/g c ) 2 


4 (g - \q\)/q c 



Table 1. The scaling exponent r(q) for the weak and strong disorder regimes. 

Computations were done through the Monte-Carlo method on a 600 x 600 lattice. The 
simulation consisted of 10 3 Monte-Carlo steps per site. The averages were taken on a 100 x 
100 sublattice. Weak and strong regimes were investigated, corresponding, respectively, to 
q c = 10 and q c = 1. Here, we report data for q = 3, 6, 9 in the weak regime, and q = —3, —6 
in the strong regime. We considered a single realization of (f>(x) for each pair (q, q c ) of 
parameters. As a matter of fact, due to the self-averaging property of model (2.1) [14], a 
single realization of the multifractal wavefunction leads to scaling exponents identical to the 
ones computed from quenched averages in the large box limit. The numerical results are 
shown in Figs. 1 and 2, toghether with the exact slopes for comparison. We find in a very 
clear way that the weak regime is well-accounted for ESS-I, while the strong regime exhibits 
ESS-II. 

III. EXTENDED SELF-SIMILARITY IN STRANGE ATTRACTORS 

There is an element in the cascade description of multifractal wavefunctions proposed 
above that brings our attention to the self-similar structure of strange attractors: the es- 
sential point is that the wavefunction normalization is preserved throughout the cascade. 
This simple constraint makes a bridge to the phenomenological analysis of strange attractors 
established by Benzi et al. [8]. In their work, a large and fixed number of N points homoge- 
neously distributed in some compact region of space (the basin of attraction) is sucessively 
iterated, producing an approximation to the attractor. It is assumed a set with the same 
statistical properties would be generated after N iterations of any arbitrary initial point, for 
N large enough. Structure functions of order q may be defined as 

q 

(3.1) 

where n(r) is the number of points contained in a ball of radius r, while d^j) is the distance 
between points labelled by i and j, and Q(x) is the Heaviside function. As N — > oo, 
multifractality means that S q (r) ~ at small scales. Since the total number of points 
N is kept constant through iterations, the quantity n(r) is immediately recognized as the 
analogous of J n d 2 x\ifj\ 2 in the context of localization. The main difference between the 
dynamics of strange attractors and the evolution of delocalized wavefunctions, however, is 
the fact that the latter will not have a finite area support in general. 

The natural question we pose regards the corrections to the multifractal asymptotic 
behavior at large scales, comparable to the overall size of the attractor. We performed 
numerical experiments on the Henon and Lorenz strange attractors. The strange attractor 
generated by the Henon mapping [16], 

Xn+l 1 ax n -\- y n , 

y n+ i= bx n , (3.2) 



S g (r) 



(n(r) q ) 



N 



N 

E 

i=i 



N 
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with a = 1.4 and b = 0.2, is fully chaotic, containing no periodic orbits. As it is well- 
known, at certain range of scales the Henon attractor is locally like several approximately 
parallel stripes. At large enough scales, we may say, therefore, Henon attractor is roughly 
one-dimensional. However, sweeping a considerable larger range of scales, the fractal sub- 
structure shows up and one gets a non-trivial fractal dimension slightly greater than unity. 
To study such a dimensional crossover, we have considered the Henon mapping for 10 5 
iterations, with initial point x = y = 0. We focused our attention on the 2428 points 
contained in the region \x — 0.76| < 0.1, \y — 0.16| < 0.01. Defining now the functions 
G q (r) = \n[S q (r)] — q In r, we clearly observe, from the results depicted in Fig. 3, for q = 3, 6 
the existence of ESS-II in the Henon attractor. The plateaus observed at larger values of r 
confirm the existence of one-dimensional large scale structures. 

Similar computations were carried out for the Lorenz attractor [17]. This dynamical 
system is given by the following coupled differential equations, 

x= - s (x - y) , 
y— —xz + rx — y , 

z— xy — bz , (3.3) 

where we take s = 10, b = 8/3 and r = 28. The discrete flow (x(t n ),y(t n ), z(t n )) which 
generates the strange attractor in the three-dimensional phase space is obtained using the 
time step 5 = 0.01 and iV = 3000. The initial point has coordinates x = y = and 
z = 0.01. The first 2000 iterations were discarded. The range of scales investigated runs 
from the smallest ones to the overall size of the attractor, so that we expect now a crossover 
between a set with non-trivial scaling behavior and a set of vanishing dimension (the Lorenz 
attractor as "viewed from the infinity"). We confirm ESS-I from Fig. 4, where the plots of 
G q {r) = \n[S q (r)\ for q = 3 and q = 8 are shown. Note that N is not a very large number, 
usually the condition to get good evaluations of the scaling exponents. The situation is 
here very similar to the one found in turbulence, when one extracts, through ESS, scaling 
exponents out of low Reynold's numbers data. 

We have material evidence, therefore, to hypothesize that ESS is a property shared by 
other strange attractors and metal-insulator transitions. It is not clear, however, under 
which conditions ESS-I or strict ESS-II will hold. Next, we attempt at a phenomenological 
description of ESS in general, under the light of random /9-models. 

IV. RANDOM /3-MODEL ANALYSIS 

The existence of ESS in multifractal systems seems to be intimately related to the dynam- 
ics of fragmentation. The random /3-model approach to strange attractors and turbulence 
[8], originally devised to capture the statistical properties associated with fragmentation 
events, is an interesting ground to get some clues on ESS. 

The multifractal set given by a strange attractor is generated, in the random /3-model 
line of thought, from the fragmentation of boxes in some <i-dimensional phase-space of a 
dynamical system. At the n-th fragmentation step, a box with volume Q n = £f v containing 
a certain number of points N n is fragmented into sub-boxes with volumes — ^( n +i)> 

containing each of them A( n+1 ) points. The total volume of the "newborn" boxes is a fraction 
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f3 n (a random variable) of the original volume of the "parent" box. In the random /5-model 
we assume that at each step of fragmentation the new length scale is a fixed fraction of the 
previous one, that is, £( n +i) — £ n / a i where a > 1. Furthermore, the volume ratio (5 n varies 
from box to box, without any correlation, being described by some probability function 
P(/3). Since the total number of points is conserved, we may write 



N (n+1 



A simple relation may be obtained from (4.1) for the number of points in a box generated 
in the sequence of fragmentations (3q — > f3± — > ... — > fl( n -i)' 



d 



n-l 



^ = (?l LU" 1 (4.2) 



i=0 



Thus, we find that 




where 



C(«) 

SM = (K) = I ? I , (4.3) 



C{q)=qd-\og a (p-«) . (4.4) 



Similar arguments may be advanced for the structure functions in turbulence, taking into 
account the energy cascade, giving ((q) = q/3 — \og a {f3^ l ^ q ^' > ) . 

Observe that for a given a > 1, there is a discrete set of a d possible values of (5 < 1. In 
order to approach metal-insulator transitions, we introduce a variant of the random /3-model 
which allows not only a larger set of /3's, but also to have (3 > 1. In this way, statistical 
fluctuations of the random variable j3 can be determined, in good approximation, from some 
probability density function p([3) with < (3 < oo. To define the alternative model, consider 
the union Uq of Mq boxes of volumes £q, containing each iVo points. We call Uq a "coherence 
region", since after the first fragmentation step, by imposition, a region Uoi is produced, the 
union of M 01 boxes of volumes if, each one taking the same number of N 01 = M N /M 01 
points. We have j3 = M m / 'M a~ d = N /N ia~ d , which may be now larger than unity. The 
model is supplemented by a "splitting rule" (deterministic or not) that states if Uqi is still 
a coherence region or if it will be broken into a number p of independent sub-coherence 
regions t/ i(i), t^oi(2), foi(p)> & U °f them having identical volumes. The same procedure is 
then inductively repeated as fragmentation proceeds. If, for instance, coherence regions are 
initially split into two sub-regions, we get the following diagram of fragmentations: 




^01(1)2(1) •• 

^01(1)2(2) •• 

^01(2)2(1) •• 

^01(2)2(2) •• 
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After n iterations, a box of volume £ d will contain a number of points that is still given by 
(4.2). Since all paths in the fragmentation diagram have equivalent probabilistic weights, 
for reasonable choices of the splitting rule, then it follows that both (4.3) and (4.4) hold as 
well, if the number of coherence regions is large enough to allow for statistical averaging. 

As an application of the above ideas, we identify N n /N with the box probability J Qn \tp\ 2 
in condensed matter localization. Recalling the discussion of sec. II, we reproduce the 
scaling exponents for the ranges \q\ ^> q c in the weak and strong regimes of disorder in 
terms of a bifractal model. Let us take /3_ < 1 < j3 + as the volume ratios, which are 
generated with probabilities P_ and P + = 1 — P_, respectively. When q — > ±oo, it folows 
that (f3~ q ) — > P T j3^ q . We find that both disorder regimes correspond to P + = P_ = 1/2 
and a = y^2. Hovewer, the weak disorder regime gives 

log 2 /3 ± = -(-±2) , (4.5) 
qc qc 

while the strong disorder regime gives 

log 2 /3_ = -l, log 2 /3+ = i-l. (4.6) 

Qc 

The phenomenon of ESS reflects, in loose words, the fact that some fundamental sta- 
tistical property of the self-similar cascade keeps holding in the crossover towards non- 
multifractal behavior, while a set of parameters (as the fragmentation parameter a or the 
dimension d) evolves "adiabatically" . Let us assume that the probability density of (3 is 
fixed along the cascade, in the random /3-model framework. We are able to establish, in this 
way, two independent phenomenological pictures that support ESS. 

Picture A: Target Space with Fractal Dimension. We assume in this case that frag- 
mentation, initially defined in d dimensions, has been gradually modified to occur in a 
(i/a-dimensional space. Considering that the probability density determining the relative 
number of boxes generated within coherence regions does not vary, the only way to preserve 
the statistics of f3 is from a modification of the parameter a, as may be inferred from (4.1). 
We have 

Pn = = ^(aT d/a (4.7) 

and, therefore, a' = a a . Performing the substitution d — > d/a and a — > a a in the expression 
(4.4), we obtain ((q) — > ((q)/a for strange attractors or in the localization problem, which 
implies ESS-I in these systems. 

Picture B: Enhanced/ Supressed Fragmentation. We consider here the situation where 
fragmentation is enhanced or supressed. That is, the relative number of boxes generated 
within coherence regions gets a factor c. The fragmentation takes place in <i-dimensional 
space. As in the previous picture, the only way to keep the probability density for (5 fixed 
is from a modification of the parameter a. We have 

Pn = = • (4.8) 

Defining a from c = a d<yl ~ a \ we get, from (4.8), a' = a a . We will have now ESS-II, with the 
crossover functions given by fi(r) = r ^ 1 " 1 /") and / 2 ( r ) = r l / a . 
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An arbitrary combination of pictures A and B is the most general situation, leading 
always to ESS-II. It is interesting to note, in passing, that both of these pictures yield 
ESS-II for the turbulence cascade. 

V. CONCLUSION 

We studied infra-red effects in an exactly solvable two-dimensional metal-insulator tran- 
sition modeled by random Dirac fermions, observing in a clear way the existence of ESS-I 
and ESS-II, in the strong and weak disorder regimes, respectively To our knowledge this 
is the first verification of ESS for a system that is not manifestly classical. A cascade de- 
scription of the multifractal probability density profile was proposed along the lines of the 
random /3-model, in a spirit similar to the usual applications performed in strange attractors 
and turbulence. Furthermore, we found that strange attractors can in fact show ESS, from 
a straightforward numerical analysis of the Henon and Lorenz dynamical systems. 

Two independent pictures, defined within the random /3-model approach, were put for- 
ward as possible phenomenological descriptions of ESS. A very interesting (and challenging) 
problem is to check, then, if these pictures are somehow realized in the multifractal systems 
where ESS holds. We believe, however, that the next natural step of research concerns the 
issue of ESS in other metal-insulator transitions already known in condensed matter. It is 
worth mentioning that in the two-dimensional localization problem, r(l) = £(l) — 2 = 
is an exact relation 2 , in perfect analogy with Kolmogorov's 4/5 law of turbulence. There- 
fore, it is likely ESS can play for the subject of localization the same important role it 
has for turbulence, providing better evaluations for the multifractal exponents of structure 
functions. 
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2 The exact result = d in ci-dimensional space follows from translation invariance of the box 
probabilities. In the language of the random /3-model, we have, according to (4.4), ((3) = 1. 
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FIGURES 



FIG. 1. 

ESS-II data for structure functions in the weak disorder regime (q c = 10), taking results 
for q = 3, 6, 9. 

FIG 2. 

ESS-I data for structure functions in the strong disorder regime (q c = 1), taking results 
for q = —3, —6. 

FIG 3. 

ESS-II data for the Henon attractor, where structure functions of order q = 3, 6 are taken 
into account. 

FIG. 4. 

ESS-I data for the Lorenz attractor, where structure functions of order q = 3, 8 are taken 
into account. 
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